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Motion of Major Ice Shelf Fronts in Antarctica from Slant Range Analysis of 
Radar Altimeter Data, 1978 - 1998 

H. J. Zwally, M. A. Beckley, A. C. Brenner, and M. B. Giovinetto 

Slant range analysis of radar altimeter data from the Seasat, Geosat, ERS-1 and ERS-2 databases 
are used to determine barrier location at particular times, and estimate barrier motion (km/yr) for 
major Antarctic ice shelves. The barrier locations, which are the seaward edges or fronts of 
floating ice shelves, advance with time as the ice flows from the grounded ice sheets and retreat 
whenever icebergs calve from the fronts. The analysis covers various multiyear intervals from 
1978 to 1998, supplemented by barrier location maps produced elsewhere for 1977 and 1986. 
Barrier motion is estimated as the ratio between mean annual ice shelf area change for a 
particular interval, and the length of the discharge periphery. This value is positive if the barrier 
location progresses seaward, or negative if the barrier location regresses (break-back). Either 
positive or negative values are lower limit estimates because the method does not detect 
relatively small area changes due to calving or surge events. The findings are discussed in the 
context of the three ice shelves that lie in large embayments (the Filchner-Ronne, Amery, and 
Ross), and marginal ice shelves characterized by relatively short distances between main 
segments of grounding line and barrier (those in the Queen Maud Land sector between 10. 1 W 
and 32.5° E, and the West and Shackleton ice shelves). All the ice shelves included in the study 
account for approximately three-fourths of the total ice shelf area of Antarctica, and discharge 
approximately two-thirds of the total grounded ice area. 
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ABSTRACT 


Slant range analysis of radar altimeter data from the Seasat, Geosat, ERS-1 and 
ERS-2 databases are used to determine barrier location at particular times, and 
estimate barrier motion (km a' 1 ) for major Antarctic ice shelves. The analysis 
covers various multiyear intervals from 1978 to 1998, supplemented by barrier 
location maps produced elsewhere for 1977 and 1986. Barrier motion is 
estimated as the ratio between mean annual ice shelf area change for a 
particular interval, and the length of the discharge periphery. This value is 


positive if the earner location progresses seaward, or negative it the barrier 
location recesses '"break-base. Either positive or negative 'values a~e fewer limrt 
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estimates because the method does not detect relatively small area changes due 
to calving or surge events. The findings are discussed in the context of the three 
ice shelves that lie in large embayments (the Filchner-Ronne, Amery, and Ross), 
and marginal ice shelves characterized by relatively short distances between 
main segments of grounding line and barrier (those in the Queen Maud Land 
sector between 010.1°W and 032.5°E, and the West and Shackleton ice 
shelves). All the ice shelves included in the study account for approximately 
three-fourths of the total ice shelf area of Antarctica, and discharge 
approximately two-thirds of the total grounded ice area. 

INTRODUCTION 


The ice shelves of Antarctica account for 1 1 % of the area and 2.5 % of the ice 
volume (Drewry, 1983a). Mass budget studies of the whole ice mass of 
Antarctica list estimates of calving that account for 80 % of total mass output with 
an assessed error that is larger than the net mass budget estimate (e.g. Orheim, 
1985; Jacobs and others, 1992, 1996). Estimates of calving normally include ice 
shelf discharge across fixed reference peripheries, and the change in location of 
the ice shelf front (barrier) of large ice shelves is one of many indicators of global 
change (e.g. SCAR / IGBP, 1 989). 

In the following sections we summarize mean annual barrier motion estimates 


otained from slant range analysis of radar altimeter data (Martin and others. 


'homas and othe-s, 1383: Zwaily and others. ' 
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Zwally, 1989; Zwally and Brenner, 2000) selected from the Seasat, Geosat, ERS- 
1 and ERS-2 databases covering diverse periods between 1978 and 1998, 
supplemented by barrier location maps produced elsewhere for 1977 and 1986. 
As described in the following section, our estimates are the lower limit either of 
seaward barrier motion or of barrier break-back (regression) for the observation 
periods because the method does not account for either calving or surge events 
smaller than the barrier motion. The ice shelves selected for our study comprise 


approximately three-fourths of the total ice shelf area of Antarctica, drain 
approximately two-thirds of the total grounded ice area, and are representative of 
the two principal types of ice shelf (e.g., Giovinetto, 1970) (Figure 1). Ice shelves 
of the first type lie in large embayments and are characterized by long distances 
between main segments of grounding line and barrier, as well as large areas of 
subglacial net melting and freezing. We include all three: the Filchner-Ronne, 
Amery and Ross ice shelves (FRIS, AIS, and RIS, respectively), each ice shelf 
split in three lobes, each lobe corresponding to a different drainage system. Ice 
shelves of the second type lie either in relatively small embayments or open, and 
are characterized by relatively short distances between main segments of the 
grounding line and barrier, as well as subglacial net melting. We include those 
that lie in the Queen Maud Land sector (QMLs), and the West and Shackleton 
ice shelves (WIS and SIS, respectively). For purposes of discussion, we created 
sectors by selection of prominent capes inland in the QMLs, and a major ice rise 


in the WIS, and projecting a line along their particular longitude 
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(the capes and ice rise are clearly depicted on Figure 1). Ice shelf lobes and 
sectors are designated western, central and eastern (W, C, and E, respectively). 


BARRIER LOCATION 

The slant range analysis of radar altimeter data to compile barrier location is 
based on the short time (of the order of one second) during which the radar 
altimeter that detects backscatter signals on a small surface-elevation range or 
window, fails to adjust to an abrupt elevation change (Martin and others, 1 983; 
Zwally and Brenner, 2001). During this time the altimeter continues to detect the 
signals from a surface at the same level (e.g. from sea or sea ice in front of a 
barrier) but at an increasing range that is slanted backward (Thomas and others, 
1983). In the reverse sequence, from ice shelf to ocean, the slant range 
measurement continues backward to the higher ice shelf elevation at the barrier 
location that is closest to the satellite. The slant range analysis produces a 
distinct “V” pattern of selected reflection points (Figure 2) that is symmetrically 
distributed relative to the ground track (Zwally and others, 1987). This pattern, 
compiled for a series of orbits, is used to delineate the barrier location by 
manually joining the one arm of the “V” that best aligns with the orientation of the 
barrier as known from other satellite imagery, or surveys made from aircraft or 
ships. The approach can not be applied to all barrier segments, particularly if 
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approximately 50 km or larger, or if total barrier length is of approximately 50 km 
or less. The databases used in this study produced compilations of “V” patterns 
per 100 km of discharge periphery ranging from 29 for the WIS, to 69 for the RIS. 
These average values varied widely for particular lobes or sectors, and 
databases, of any one ice shelf; some examples of poor coverage are mentioned 
in following sections. 

The procedure to delineate a barrier is relatively simple in sectors where barrier 
alignment is characteristically smooth as observed for long periods (e.g. sectors 
of the Filchner-Ronne and Ross ice shelves). The procedure becomes more 
difficult if barrier alignment is either characteristically uneven, or varies in 
relatively short periods, or both (e.g. sectors of the West and Shackleton ice 
shelves). Where it is not possible to rely on the alignment of one arm consistent 
with the known orientation of a barrier, the barrier location is inferred by joining 
the vertices of “V” patterns. In these cases there is a larger probability of 
erroneously including in the delineation one or more patterns produced by an 
iceberg lying close to the barrier, or ice shelf features that are not the barrier (e.g. 
rifts and ice rises; Stephenson and Zwally, 1989). Once completed, a particular 
barrier delineation is assigned the center date of the period between the first and 
last orbit used to compile “V” patterns. 


In general, for each orbital path and altimeter waveform data set used in the 
compilation, there are several other sets that were examinee ana oiscameo 
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because the degraded waveforms were not suitable for analysis. A degraded 
waveform set is generally associated with a more oblique angle between orbital 
path and barrier alignment. Missing data results in longer distances between 
plotted “V” patterns and thus greater interpolation distance (the interpolation 
procedure is illustrated in Figure 2). However, a greater reliance on waveform 
sets associated with orbital paths in which the ground track approaches the 
barrier at or near a right angle, significantly reduces the error in the determination 
of a barrier location. The error in the determination of a particular “V” pattern 
location decreases from + 1 km to + 0.1 km as the cross over angle is close to 
perpendicular (Thomas and others, 1983; Zwally and others, 1987; Stephenson 
and Zwally, 1989). 

BARRIER MOTION 


Barrier motion (M, in km a' 1 ) is estimated as the ratio between ice shelf area 
change (D, in km 2 ) determined for a particular interval (T, in a), and the length of 
the discharge periphery (L, in km): 

M = (D T' 1 ) L' 1 (1) 


In Figure 3 it is shown that the length of the discharge periphery is measured 

along a smooth line traced between the two delineated barrier locations. If D is a 
positive number (rtgure 3. A), M (> 0) is the lower limit of barrier motion seaward 
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because calving events smaller than the area change in particular sectors are not 
detected. If D is a negative number (Figure 3.B), M (< 0) is the lower limit of 
barrier regression (break-back) because seaward barrier motion smaller than the 
area change in particular sectors (e.g. due to prior surges of outlet glaciers or ice 
streams) are not detected. 

Any fortuitous detection of calving (Figure 3.A) where the delineation of the 
barrier for the interval closing date lies land ward of the delineation for the 
interval onset date, generates a partial area of negative value (-d, in km 2 ) that is 
factored by -1 and added to D: 

M > 0 = ({D + [(-d) (-1)]} T’ 1 ) L' 1 (2) 

The estimate of M remains a lower limit estimate because -d is a measure of 
minimum detected calving. Similarly, any fortuitous detection of seaward barrier 
motion greater than the break-back rate (Figure 3.B) results in the delineation of 
the barrier for the interval closing date to lie seaward of the delineation for the 
interval onset date, generating a partial area of positive value (d, in km 2 ) that is 
factored by -1 and added to -D: 

M < 0 = {{-D + [(d) (-1 )]) T' 1 ) L' 1 


( 3 ) 
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The estimate of M remains a lower limit estimate because d is a measure of 
minimum detected seaward motion. 


DATA COMPILATION AND ANALYSIS 


Selection of ice shelves, ice shelf sectors, and barrier segments 


The data on barrier motion corresponding to the FRIS, AIS, and RIS are 
compiled for lobes that correspond to separate drainage systems (Table 1 , 

Figure 4). The criteria for the selection of drainage divides will be discussed in 
detail elsewhere; however there are several aspects that provide a perspective to 
the findings listed in following sections. Principally, it should be noted that the 
new divides pattern is only broadly similar to that discussed in prior studies 
(Giovinetto and Bentley, 1985; Bentley and Giovinetto, 1991; Vaughan and 
others, 1999; Giovinetto and Zwally, 2000). The new pattern is based on the 
ERS-1 radar altimeter database transformed in a 5 km grid format listing 
maximum surface slope gradient and orientation (Zwally and Brenner, 2001). 
Southward of 81.5°S (outside ERS-1 and ERS-2 radar altimeter coverage) the 
grid was extended using a digital elevation model (Bamber and Huybrechts, 

1996) based on the Scott Polar Research Institute map (Drewry, 1983b). The 
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separated from the ice shelf flow that originates in part as land ice in East 
Antarctica (systems 2, 3 and 17). 

Splitting the ice flow of the FRIS and RIS in the context of land ice provenance 
requires further analysis in the delineation of divides southward of ERS-1 and 
ERS-2 radar altimeter coverage, particularly segments of divides on the 
grounded ice sheet between systems 1-2-3, and 17-18-19. Other aspects and 
features of the divides pattern require more study. For example, ice flow from the 
area of subglacial Lake Vostok (centered at approximately 77°S, 105°E; cf. 
Zwally and Brenner, 2001) is shown in Figure 1 as part of drainage system 17, 
although a flow lines pattern derived from the Radarsat database (K. Jezek, 
personal communication, 2001) shows it as part of system 13. The delineation of 
divides were extended on ice shelves following ice flow boundaries selected from 
compilations for the FRIS (Swithinbank and others, 1987), AIS (Hambrey and 
Dowdeswell, 1994), and RIS (Thomas and others, 1984; Fahnestock and others, 
2000). The new divides do not coincide with the ridges of Berkner Island (Bl, in 
system 2) and of Roosevelt Island (Rl, in system 19), and the coastal end point 
of the divide between systems 16 and 17 has been shifted from Cape Crozier 
(CC) in Ross Island, to the northernmost point on Brown Peninsula (77.78°S, 
165.6°E). 


It should also be noted that because radar altim 

narrow coastal zone characterized by complex to 
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gradient, the divides were extended to the coastline of the grounded ice sheet 
based on surface topography maps (Drewry, 1983b; Ferrigno and others, 1996). 
Most grounding lines and junction points (where the coastline, grounding line and 
barrier adjoin) used as guides for some divides were those compiled in a 
separate map (Swithinbank, 1988). We introduced a few changes due to large 
scale calving (e.g. the divide between drainage systems 12 and 13 was shifted 
westward by approximately three degrees of longitude due to large scale calving 
in the SIS). Lastly, system 10 now includes the three major glacier basins in the 
interior, i.e. the basin area of the Fisher Glacier, formerly part of system 9, is 
joined to the Mellor and Lambert glacier basins. 

Filchner-Ronne Ice Shelf 

The bulk of the slant range data correspond to the FRIS-W and -E lobes using 
ERS-1 and ERS-2 databases covering an interval of 5.91 a between 1992-06-08 
and 1998-05-05. The FRIS-C lobe presents a special case because ice flow is 
split by Berkner Island, and it was necessary to consider three separate barrier 
segments designated Cl , C2 and C3, from west to east, respectively. For the 
FRIS-C1 segment we used the delineation based on the ERS-1 database 
selecting orbital sets from its Geodetic Mission (GM) phase centered on 1994-09- 
26. We could not produce any other reliable delineation based on slant range 
analysis, and therefore compared it with the barrier location shown in s. map from 
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Multispectral Scanner database centered on 1986-02-21 supplemented by 
information from the EOSAT database (Swithinbank and others, 1987). We were 
able to produce delineations for the central and eastern segments (FRIS-C2 and 
-C3) using the ERS-1 and ERS-I(GM) databases centered on 1992-06-08 and 
1 994-09-26, respectively. 

The findings for the FRIS indicate a mean motion of 1 .24 km a' 1 along a 
discharge periphery of 891 km, the mean for each lobe showing small differences 
from the overall mean, i.e. the M values for the FRIS-W, -C, and -E are 1 .25, 
1.15, and 1.38 km a' 1 . The largest difference in M values is found between 
barrier segments of distinct physiographic characteristics within FRIS-C. In this 
lobe with a discharge periphery of 282 km, the M value for segment Cl is 2.13 
km a" 1 along 96 km or 34 % of the total periphery, and for segment C2 is 0.17 km 
a' 1 along 1 14 km or 40 % of the total periphery. A difference by a factor of 12.5 is 
not surprising as segment Cl is the discharge periphery for flow that includes the 
input from at least the western and central parts of the Foundation Ice Stream 
basin (assuming that flow from the eastern part of the basin might flow eastward 
of Berkner Island as suggested by detailed surface topography (cf. Sievers and 
others, 1993)), while segment C2 is the drainage periphery for the northern 
sector of the ice cap on Berkner Island. The M value for segment C3 is 1 .40 km 
a' 1 along 72 km or 26 % of the total periphery for the lobe. 
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As an example of the caution to be exercised in using the data for the FRIS and 
all other ice shelves included in this study, a major calving event (approximately 
8200 km 2 ) was recorded in 2000-05 by NOAA-14 along the westernmost 340 km 
of the FRIS-W periphery, with the maximum distance between the former and 
new barrier (break line) being slightly over 40 km in some sectors. 

Ice shelves of the Queen Maud Land sector 

We produced a single barrier delineation for the QMLs using the Geosat (GM) 
database centered on 1986-01-01 , and compared it with barrier locations shown 
in maps printed by the United States Geological Survey (USGS; Ferrigno and 
others, 1996) for sub-sectors QMLs-W and -E, and by the Scott Polar Research 
Institute (SPRI; Drewry, 1983b) for QMLs-C. The USGS compilation for QMLs-W 
is based mainly on NOAA-7 imagery obtained on 1 983-02-06, supplemented by 
imagery from NOAA-12 obtained on 1983-01-29 for a small part of the 
westernmost sector (between 006. 1°W - 003. 5°W). We dated the barrier mapping 
for the whole sub-sector at 1983-02-06 and assigned an interval of 2.90 a. The 
USGS compilation for QMLs-E is based on NOAA-7 imagery obtained on 1983- 
01-29, and we assigned an interval of 2.92 a. The SPRI compilation for QMLs-C 
is based on a report from ship and aircraft surveys completed in the summer of 
1976-1977 (Cooper and others, 1983), and we assigned an interval of 9.00 a. 
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The findings for the QMLs indicate a mean motion of -1.15 km a' 1 along a 
discharge periphery of 2415 km, an overall mean matched by the M value for 
QMLs-W (-1 .02 km a' 1 along a periphery of 1071 km) but not in sub-sectors 
QMLs-C and -E (the M values are -0.59 km a' 1 over a periphery of 964 km, and - 
2.98 km a' 1 over a periphery of 380 km, respectively). In part, the larger break 
back rate in QMLs-E might be explained by the larger land ice discharge 
expected from the eastern area of system 6 which extends farther inland than the 
western area of the system because the difference in mean net accumulation at 
the surface is small (e.g. Giovinetto and Zwally, 2000). Also in part, the smaller 
break back rate in the QMLs-C might be explained by impeded land ice 
discharge due to the mountain range alignment inland and parallel to the 
grounding line, as well as the absence of subglacial troughs, outlet ice streams, 
etc. (e.g. van Autenboer and Decleir, 1978; cf. Bentley and Giovinetto, 1991). 


Amery Ice Shelf 


The best slant range measurements set compiled for the AIS barrier was 
obtained from the Geosat(GM) database centered on 1986-01-01, and adequate 
compilations were obtained from the Seasat database centered on 1978-09-15 
and the ERS-I(GM) database centered on 1994-09-26. The M values obtained 
for all three lobes were uniform between the two possible intervals with no 


extraordinary cross over of the barrier traces, i.e. any measured d and -d values 


were of iess than +10 km 2 . Therefore, for simpibitv. we present a sincie 


assessment for the longest possible interval (16.04 a). 
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The mean motion for the AIS is estimated at 0.60 km a' 1 along a discharge 
periphery of 294 km. There is a small difference between the M values for the 
AIS-W and AIS-E lobes (0.41 along a periphery of 145 km or 49 % of the total, 
and 0.55 km a' 1 along a periphery of 80 km or 27 % of the total, respectively). 
These two lobes drain the grounded ice areas extending on either side of the AIS 
(systems 9 and 1 1 ). As expected, the largest M value is estimated for AIS-C 
(1 .03 km a' 1 along 69 km or 24 % of the total periphery) that drains the three 
large interior basins of the Fisher, Mellor and Lambert glaciers (system 10). 

West and Shackleton ice shelves 

The barrier delineation for both the WIS and SIS was based on the Geosat(GM) 
database centered on 1986-01-01 . These barrier locations were compared with 
those obtained from the USGS map, first for the WIS-W sector based on NOAA-6 
data centered on 1980-01-18 (an interval of 5.95 a), and second for the WIS-E 
sector and SIS, both based on NOAA-7 data centered on 1983-02-16 (an interval 
of 2.87 a). 

The mean motion for the WIS is estimated at -1 .19 km a' 1 along a discharge 
periphery of 605 km. The two component estimates for the WIS-W and -E 

sectors show a relatively large range of M values given the relatively uniform 
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Smith, 1985; Budd and Warner, 1996), suggesting that our sampling intervals of 
different length (-0.93 km a' 1 for 5.95 a, and -1.49 km a' 1 for 2.87 a, respectively) 
might include or exclude calving or surge phenomena on either sector that 
occurred between January 1980 and January 1986. 

Ross Ice Shelf 

Slant range measurements from the ERS-1 and ERS-2 databases centered on 
1 992-06-25 and 1 998-05-05, respectively, produced barrier locations 
encompassing an interval of 5.86 a that appears to be free of major calving 
events or surge episodes. The boundaries between lobes, inland as well as on 
the ice shelf, correspond to boundaries between ice streams typically identified at 
their gates near the grounding line of western West Antarctica. It was mentioned 
above that the coastal end point of the divide between systems 16 and 17 is the 
northernmost point on Brown Peninsula (77.78°S, 165.6°E). In this study, for 
brevity, we do not include the minor RIS flow southward of Ross Island into the 
area of the McMurdo Ice Shelf (cf. Giovinetto and Zumberge, 1968; Swithinbank, 
1970), and select Cape Crozier on Ross Island as the western “boundary” for the 
RIS-W lobe. Elsewhere we use the divides as boundaries; the boundary between 
systems 17 and 18 is the boundary between ice streams A and B, and that 
between systems 18 and 19 is the boundary between ice streams C and D (cf. 
Thomas and others, 1984; Fahnestock and others, 2000). 
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Relative to the mean motion for the RIS as a whole (1 .37 along a discharge 
periphery of 861 km), the larger M value estimates are for the RIS-W and -C 
lobes (1.44 km a' 1 along 327 km or 38 % of the periphery, and 1.84 km a' 1 along 
243 km or 28 % of the periphery, respectively). The lower M value estimate is for 
the RIS-E lobe (0.91 km a' 1 along 291 km or 34 % the periphery). The discharge 
from this lobe is split by Roosevelt Island in complex interaction with the ice flow 
from ice streams D and E (Fahnestock and others, 2000). 

DISCUSSION AND CONCLUSIONS 

Our determination of motion based on area change overtime produces lower 
limit estimates because the approach does not detect area change due to 
relatively small calving or surge events. This shortcoming is shared with all other 
methods that estimate barrier motion based on differences in barrier location 
overtime (e.g. Keys and others, 1998). However, the area change approach is 
not affected by the inherent sources of error in estimates of ice shelf motion 
based on vector analysis, namely velocity interpolation between survey sites that 
are far apart, and azimuth variation of the motion (e.g. Hofmann and others, 

1964; Lisignoli, 1964). Preceding studies show large directional (lateral) 
variability of ice shelf motion for periods of a few years and longer, particularly in 
the proximity of either permanently or temporarily grounding areas (e.g. Lisignoli, 
1964; Giovinetto and Zumberge, 1970; Swithinbank and others, 1987; Keys and 
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We assess the composite error for each estimate of motion produced in this 
study on the basis of the mid-range value for the error in the compilation of 
barrier location from slant range analysis (±0.5 km), and of the mean error in the 
location of features, including barriers, applicable to the supplementary maps 
used in the study (±2.5 km; e.g. Kohnen, 1982; Ferrigno and others, 1996). 
Treated as standard errors, the composite values are ±0.71 km for estimates 
based on two barrier locations determined from slant range analysis, and ±2.55 
km for estimates based on one of two barrier locations determined from maps 
compiled elsewhere. These two composite values are allocated to the 
corresponding estimates of motion listed in Table 1, split over each particular 
interval. 


Excluding the extraordinarily large error applicable to the estimate for FRIS-C2, 
the assessed error ranges from ±4 % for AIS-C (1 .03 ± 0.04 km a' 1 ), to ±86 % for 
QMLs-W (-1 .02 ± 0.88 km a' 1 ). The error is generally smaller for the lobes 
corresponding to ice shelves in large embayments, where it is between the 
aforementioned ±4 % for AIS-C, and ±26 % for FRIS-C (-1.15 ± 0.30 km a ’), 
than for the sectors in marginal ice shelves, where it is between ±29 % for QMLs- 
E (-2.98 ± 0.87 km a 1 ), and the aforementioned ±86 % for QMLs-W. Both the 
slant range method used by us and the imagery utilized by others to compile the 
supplementary maps produce excellent spatial coverage relative to vector 


>!\'sis, and significantly reduce the probable s-ror in the determir.atio: i of 
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motion. Nonetheless, the errors listed in Table 1 are large and suggest that 
errors listed in other studies may be grossly underestimated. 

Our findings are not directly comparable to many estimates of motion reported in 
the literature because the majority are based on vector analysis, as well as in 
surveys normally made at distances ranging from a few to more than 100 km 
from the barrier. Nevertheless, the following comparisons provide a perspective 
on the results listed in Table 1 : 

(a) Motion for a sector equivalent to FRIS-C3 and FRIS-E combined has been 
estimated from three vectors for sites a few kilometers from the barrier at 1 .33 
km a 1 for a period ca. 1957-1962 (Lisignoli, 1964). Our weighted estimate of 
motion is 1 .39 + 0. 1 8 km a' 1 (ca. 1 992-1 998). 

(b) Motion for a sector equivalent to the central and eastern parts of AIS-C has 
been estimated at 0.98 km a' 1 at a mean distance from the barrier of 
approximately 37 km for (the range for the velocity and distance of the nine 
vectors ca. 1 962-1965 is from 0.81 km a' 1 at 59 km from the barrier to 1 .25 km a 
1 at 15 km from the barrier; Budd and others, 1982). This estimate suggests a 
mean motion of approximately 1 .3 km a* 1 at the barrier, a rate confirmed for the 
central part of the AIS for 1 997 (Fricker and others, this issue). Our estimate for 
the lobe, including its western part, is 1 .03 + 0.04 km a' 1 (ca. 1978-1994). 

(c) Mean motion for the equivalent of RIS-C can be estimated at 0.84 km a' 1 for 
ca. 1962-1966 from 28 vectors evenly distributed along a discharge periphery of 


f 
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(Giovinetto and Zumberge, 1968; original data from a personal communication by 
E. Dorrer, 1966 (cf. Hofmann and others, 1964; Dorrer, 1970)). This estimate 
suggests a mean motion of approximately 1 .0 km a' 1 at the barrier, which is well 
under our estimate of 1 .84 + 0.12 km a' 1 (ca. 1992-1996) for a discharge 
periphery of exactly the same length (see ‘e’ below). 


More direct comparisons are possible with information obtained from reports that 
either show different barrier location for particular periods, or produce estimates 
from the same type of compilations: 

(d) Detailed mapping of the barrier location in the FRIS (Swithinbank and others, 
1987) shows the barrier location for the equivalent of FRIS-W during February- 
May 1986, and its location during February 1973 for the western 1/3 of its length, 
and February 1974 for the other 2/3 of its length. The mean distance between the 
two barrier traces is approximately 14 km, from which we obtain a weighted 
mean motion of 1 . 1 4 km a' 1 . Our estimate of motion for the same lobe is 1 .25 + 
0.10 km a' 1 (ca. 1992-1998). 

(e) A compilation of the RIS barrier location for particular years distributed 
between 1841 and 1997 includes reports on detailed mapping from ships (GPS 
navigation and radar detection of the barrier) for the years 1987, 1994, and 1997, 
when the location was determined with an error of +0.25 km (Keys and others, 
1998). Coincidentally during those years, the sector equivalent to our RIS-C lobe 


was least affected by calving or surge events (as opposed to phenomena in the 
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periods in the sector equivalent to the RIS-C lobe is relatively steady at 
approximately 0.9 km' 1 (Keys and others, 1998). Our estimate of motion for the 
lobe is 1 .84 + 0. 1 2 km a' 1 (ca. 1 992-1 998). At the time of writing we do not have 
an explanation for the large difference between the two estimates; the slant 
range data sets are among the best in distribution and quality, there is no 
evidence of large calving or surge events, and we are not aware of any reports 
on transverse rift formation. A rate of motion similar to our estimate for the RIS-C 
lobe has been reported for the central part of the RIS-W lobe and obtained by the 
procedure based on change in barrier location ca. 1983-1987 (approximately 1.9 
km a* 1 , citing the widening of existing transverse rifts and formation of an 
additional rift during the period; Keys and others, 1998). However, for the same 
period and using the same procedure, the motion for the equivalent to our RIS-C 
is approximately 0.6 km a' 1 (Keys and others, 1998), i.e. a 33 % lower value than 
the estimate of 0.9 km a-1 for 1992-1998. The estimates of motion based on 
vector analysis support the estimate of approximately 0.9 km a' 1 rather than our 
estimate of 1 .8 km a' 1 . Firstly, there is our own estimate of motion from the 28 
sites surveyed ca. 1 962-1 966 (approximately 1 .0 km a' 1 at the barrier) as 
described in ‘c’, above. Secondly, we estimate a velocity of 0.96 km a' 1 (ca. 
1976-1977) for the sector from a set of five survey sites located from a few to 35 
km from the barrier (Thomas and others, 1984; sites Q9, and RIO-13) that also 
suggests a barrier motion of approximately 1 .0 km a' 1 . More importantly, our own 
database shows an even larger difference between the estimates of motion 
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1998/05/05, corresponding to slant range analysis of the ERS-1, ERS-I(GM), 
and ERS-2 databases, respectively. The first interval (2.25 a) produces an 
extraordinarily large value of 3.37 km a' 1 , while the second (3.61 a) produces a 
value of 0.89 km a' 1 that is in close agreement with the estimates of Keys and 
others (1998). At the time of writing, we do not have an explanation for the large 
discrepancy in our data either, and therefore suggest that the estimate of motion 
listed in Table 1 for the RIS-C lobe should be used with caution and in full 
awareness of the diverse estimates obtainable for the lobe. 

Our study suggests that ice shelves occupying large embayments (the FRIS, 

AIS, and RIS) may, for periods of the order of 10 years, sustain motion seaward 
(positive M values), leading to large scale calving events. It also suggests that 
marginal ice shelves characterized by relatively short distances between principal 
segments of their grounding line and barrier (i.e. practically all other Antarctic ice 
shelves) break-back more gradually and with higher frequency (negative M 
values), and therefore large calving events are less frequent. One aspect of ice 
shelf - ocean interactions might contribute to the signed motion difference just 
described. Briefly stated, it should be expected that the FRIS, AIS and RIS, due 
to their greater thickness near the barrier, break up in response to the cumulative 
effects of tide-induced flexing (wave periods of the order 10 5 sec), while the 
thinner, marginal ice shelves break up in response to the cumulative effects of 
swell-induced flexing (wave periods of the order of 10 1 sec), even if a particular 
tide or storm surge trigger the actual break. 
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Figure captions 

Figure 1 . Barrier motion is estimated for the ice shelves that occupy large 
embayments (Filchner-Ronne, Amery, and Ross ice shelves; FRIS, AIS, and 
RIS, respectively), and some of the marginal ice shelves (in the Queen Maud 
Land sector, and the West and Shackleton ice shelves; QMLs, WIS and SIS, 
respectively). The ice drainage divides (solid red lines) show the accumulation 
areas for the ice shelves (those for the WIS and SIS can be readily inferred from 
the layout of the divides in system 12). Other features: Roosevelt Island (Rl), 
Berkner Island (Bl), Cape Crozier (CC) on Ross Island. 


Figure 2. The common “V” pattern selected reflecting points (center of each 
diamond) obtained from slant range analysis of radar altimeter data is illustrated 


for a sector of the Ross 
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the ERS-1 and ERS-2 radar altimeter databases are shown relative to the 
location compiled from NOAA-6 imagery (Ferrigno and others, 1996). The dates 
of the radar altimeter data are part of the compilation periods centered on 1 992- 
06-25 (ERS-1) and 1998-05-05 (ERS-2) mentioned in the text. 

Figure 3. Estimate of barrier motion based on the difference in ice shelf area (D) 
for an interval defined by onset and closing dates (Equation 1). The discharge 
periphery is delineated as a smooth trace between barrier locations at the onset 
and closing dates for the interval (its length = L). As described in text, the 
estimates of barrier motion do not account for all possible calving and surge 
events. (3.A) Minimum rate of barrier motion seaward, where “-d” indicates the 
minimum detected area of calving (Equation 2). (3.B) Minimum break-back rate 
(barrier regresses land ward), where “d” indicates the minimum detected barrier 
motion seaward (Equation 3). 

Figure 4. Schematic representation of barrier motion estimates listed in Table 1 . 
Solid bars indicate minimum mean barrier motion seaward, and clear bars 
indicate minimum mean break-back (barrier location regresses landward). 
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mi mm Barrier location at interval onset date 
■ ■ ■ ■ ■ Discharge periphery (trace used to measure length: L) 
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Table 1 . Summary of minimum barrier motion estimates 


Ice shelf entity 

Sector 

(degrees) 

Database 

Interval 
Onset - Closing 

(a) 

Area change i 
Diff. Partial 

f km 2 ) 
Total 

Periphery 

(km) 

Motion 
(km a' 1 ) 

FRIS-W 

061.5W-048.9W 

ERS1 - ERS2 

920608-980505 

5.91 

2977 

311 

3288 

445 

1.25 ±0.12 

FRIS-C 3 

048.9W-041.5W 

(see below) 3 


1.00* 

324* 

— 

324* 

282 

1.15 ±0.30 

FRIS-E 

041.5W-035.5W 

ERS1 - ERS2 

920608-980505 

5.91 

611 

726 

1337 

164 

1.38 ±0.12 

FRIS-whole 

061.5W-035.5W 



— 

— 

— 

— 

891 

1.24 

FRIS-C1 

048.9W-047.1W 

IfAG - ERS1GM 

860221-940926 

8.60 

1758 

— 

1758 

96 

2.13 ±0.30 

FRIS-C2 

047.1W-043.8W 

ERS1 - ERS1GM 

920608-940926 

2.30 

44 

— 

44 

114 

0.17 ±0.31 

FRIS-C3 

043.8W-041.5W 

ERS1 -ERS1GM 

920608-940926 

2.30 

232 

— 

232 

72 

1.40 ±0.31 

QMLs-W 

009.9W-009.2E 

USGS - GeoSatGM 

830206-860101 

2.90 

-2471 

-690 

-3161 

1071 

-1.02 ±0.88 

QMLs-C 

009.2E-026.5E 

SPRI - GeoSatGM 

770101-860101 

9.00 

-4414 

-659 

-5073 

964 

-0.59 ±0.28 

QMLs-E 

026.5E-032.5E 

USGS - GeoSatGM 

830129-860101 

2.92 

-3304 

— 

-3304 

380 

-2.98 ±0.87 

QMLs-whole 

009.9W-032.5E 



— 

— 

— 

... 

2415 

-1.15 

AIS-W 

070.1 E-072.6E 

SeaSat- ERS1GM 

780915-940926 

16.04 

948 

7 

955 

145 

0.41 ±0.04 

AIS-C 

072.6E-073.5E 

SeaSat- ERS1GM 

780915-940926 

16.04 

1142 

... 

1142 

69 

1.03 ±0.04 

AIS-E 

073.5E-074.9E 

SeaSat - ERS1GM 

780915-940926 

16.04 

706 

... 

706 

80 

0.55 ±0.04 

AlS-whole 

070.1 E-074.9E 



— 

... 

— 

— 

294 

0.60 

WIS-W 

081 .3E-085.6E 

USGS - GeoSatGM 

800118-860101 

5.95 

-1399 

-359 

-1758 

319 

-0.93 ±0.43 

WIS-E 

085.6E-089.6E 

USGS - GeoSatGM 

830216-860101 

2.87 

-1041 

-180 

-1221 

286 

-1.49 ±0.89 

WIS-whole 

081.3E-089.6E 



— 

— 

— 

... 

605 

-1.19 

SIS 

094.8E-102.2E 

USGS - GeoSatGM 

830216-860101 

2.87 

-5735 

-190 

-5925 

650 

-3.18 ±0.89 

RIS-W 

1 69.0E-1 79. 5W 

ERS1 - ERS2 

920625-980505 

5.86 

2758 

— 

2758 

327 

1.44 ±0.12 

RIS-C 

179.5W-169.5W 

ERS1 - ERS2 

920625-980505 

5.86 

2623 

... 

2623 

243 

1.84 ±0.12 

RIS-E 

169.5W-159.0W 

ERS1 - ERS2 

920625-980505 

5.86 

1543 

... 

1543 

291 

0.90 ±0.12 

RIS-whole 

169.0E-1 59. 0W 



— 

— 

— 

... 

861 

1.37 


a The estimate of motion for the FRIS-C lobe is the weighted mean of segments FRIS-C1, -C2, 
and -C3. 

b Large differences in estimates of motion for the RIS-C lobe are described in the text. 




